home *** CD-ROM | disk | FTP | other *** search
/ NASA Climatology Interdisciplinary Data Collection / NASA Climatology Interdisciplinary Data Collection - Disc 4.iso / software / grads / lib / cbar.gs < prev    next >
Encoding:
Text File  |  1998-04-23  |  2.0 KB  |  95 lines
  1. *
  2. *  Script to plot a colorbar
  3. *
  4. *  The script will assume a colorbar is wanted even if there is 
  5. *  not room -- it will plot on the side or the bottom if there is
  6. *  room in either place, otherwise it will plot along the bottom and
  7. *  overlay labels there if any.  This can be dealt with via 
  8. *  the 'set parea' command.  In version 2 the default parea will
  9. *  be changed, but we want to guarantee upward compatibility in
  10. *  sub-releases.
  11. *
  12. function colorbar (args)
  13. *
  14. *  Check shading information
  15. *
  16.   'query shades'
  17.   shdinfo = result
  18.   if (subwrd(shdinfo,1)='None') 
  19.     say 'Cannot plot color bar: No shading information'
  20.     return
  21.   endif
  22. * 
  23. *  Get plot size info
  24. *
  25.   'query gxinfo'
  26.   rec2 = sublin(result,2)
  27.   rec3 = sublin(result,3)
  28.   rec4 = sublin(result,4)
  29.   xsiz = subwrd(rec2,4)
  30.   ysiz = subwrd(rec2,6)
  31.   ylo = subwrd(rec4,4)
  32.   xhi = subwrd(rec3,6)
  33.   xd = xsiz - xhi
  34. *
  35. *  Decide if horizontal or vertical color bar
  36. *  and set up constants.
  37. *
  38.   if (ylo<0.6 & xd<1.0) 
  39.     say "Not enough room in plot for a colorbar"
  40.     return
  41.   endif
  42.   cnum = subwrd(shdinfo,5)
  43.   if (ylo<0.6 | xd>1.5)
  44.     xl = xhi + xd/2 - 0.4
  45.     xr = xl + 0.2
  46.     xwid = 0.2
  47.     ywid = 0.5
  48.     if (ywid*cnum > ysiz*0.8) 
  49.       ywid = ysiz*0.8/cnum
  50.     endif
  51.     ymid = ysiz/2
  52.     yb = ymid - ywid*cnum/2
  53.     'set string 1 l 5'
  54.     vert = 1
  55.   else
  56.     ymid = ylo/2
  57.     yt = ymid + 0.2
  58.     yb = ymid
  59.     xmid = xsiz/2
  60.     xwid = 0.8
  61.     if (xwid*cnum > xsiz*0.8)
  62.       xwid = xsiz*0.8/cnum
  63.     endif
  64.     xl = xmid - xwid*cnum/2
  65.     'set string 1 tc 5'
  66.     vert = 0
  67.   endif
  68. *
  69. *  Plot colorbar
  70. *
  71.   'set strsiz 0.12 0.13'
  72.   num = 0
  73.   while (num<cnum) 
  74.     rec = sublin(shdinfo,num+2)
  75.     col = subwrd(rec,1)
  76.     hi = subwrd(rec,3)
  77.     'set line 'col
  78.     if (vert) 
  79.       yt = yb + ywid
  80.     else 
  81.       xr = xl + xwid
  82.     endif
  83.     'draw recf 'xl' 'yb' 'xr' 'yt
  84.     if (num<cnum-1)
  85.       if (vert) 
  86.         'draw string '%(xr+0.05)%' 'yt' 'hi
  87.       else
  88.         'draw string 'xr' '%(yb-0.05)%' 'hi
  89.       endif
  90.     endif
  91.     num = num + 1
  92.     if (vert); yb = yt;
  93.     else; xl = xr; endif;
  94.   endwhile
  95.